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ABSTRACT 

Using two-dimensional magnetohydrodynamics calculations, we investigate a twist gen- 
eration mechanism on a magnetic flux tube at the base of the solar convection zone 
based on the idea of Choudhuri, 2003, Sol. Phys., 215, 31 in which a toroidal mag- 
netic field is wrapped by a surrounding mean poloidal field. During generation of the 
twist, the flux tube follows four phases. (1) It quickly splits into two parts with vortex 
motions rolling up the poloidal magnetic field. (2) Owing to the physical mechanism 
similar to that of the magneto-rotational instability, the roUed-up poloidal field is bent 
and amplified. (3) The magnetic tension of the disturbed poloidal magnetic field re- 
duces the vorticity, and the lifting force caused by vortical motion decreases. (4) The 
flux tube gets twisted and begins to rise again without splitting. Investigation of these 
processes is significant because it shows that a flux tube without any initial twist can 
rise to the surface in relatively weak poloidal fields. 
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1. Introduction 

Surface sunspots are the most prominent feature of solar activity. It is thought that the 
magnetic flux is gene rated by the dynamo action of the di f ferential rotation at the base of 



1955 



Choudhuri et al 



1995 



Dikpati fc Charbonneau , 



1999; 



the convection zone (iParker 

Hotta &: Yokovama , 2010alb ). The magnetic flux rises to the surface owing to the magnetic 



buoyancy and generates sunspots on the surface. 

In observational and theoretical assessments, this magnetic flux is thought to be twisted. 
On the surface, vector magnetic field obser vations have revealed t hat the sola r active regions 



have a statistically significant mean twist (jPevtsov et al 



1995 



2001 



20031 ). The quantity 



a = {Jz/ Bz) is measured, where Jz, B^, and the parenthesis denote vertical electric current, 
vertical magnetic field, and spatial average. Although the measured a has large scatter, the 
quantity has a latitu dinal dependence of a = — 2.7x lO^^^f^dcg m~^, where 0dcg is latitude in 



degree (see Fig. 2 of Lonecope et al.l . ll998r). In additio n, the theoretical (numerical) study 
of a rise in a flux tube requires twist . Ischiisslei (1979) found that an untwisted magnetic 
fiux tube quickly splits into a pair of vortex tubes of opposite circulation, which moves apart 
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horizontally and ceases to rise. This vortex is generated by t he buoyancy gradient across 
the flux tube cross-section. iMoreno-Insertis fc Emoneti (|l996f ) reported that the sufficient 



twist of a magnetic flu x tube can s uppress the splitting motion while maintaining coherency 
(see also the review of [ 



Fan 



20091) . 



There are several possibilities for generating the twist. One is that the dynamo pro- 
cess itself is responsible for it. The second possibility is that the turbulent effect stcm- 
ming from the rotation on the flux tube generates the twist, which is called the S-effect 
(ILongcope et al. ^ 19981) . Another interesting explanation for the origin of the twist is given 



by 



Choudhuri 



20031) . In his explanation, the rising toroidal field is wrapp ed by a surround 



Choudhuri et al. 



ing me an poloidal field and then gets twisted (see Fig. 4 of his 2003 paper). 
(|2004r) report that, based on this idea and the Babcock-Leighton flux-transport dynamo 
model, the latitudi nal dependence o f twist can be reproduced during most of the duration 



of the solar cycle. IChatteriee et al.l ()2006l) modeled the evolution of the rising magnetic 
flux tube with the accretion of the mean poloidal magnetic held by solving the induction 
equation. They report that without turbulent diffusivity the accreted poloidal flux is con- 
flned in a sheath at the outer periphery of the rising tube. Using plausible assumption and 
free parameters, they could obtain the value of a that is comparable to the observational 
result. 

In this paper, we investigate the generation mechanism of twist based on the idea of 
Choudhuril (|2003l ) in two-dimensional magnetohydrodynamics (MHD) calculations using 



the parameters at the base of the solar convection zone. 



2. Model 

We solve the compressive MHD equations in two-dimensional Cartesian geometry {x,z), 
where x and z denote the horizontal and vertical directions. The magnetic field has three 
components, i.e., B — {B^, By, B^) even in two dimensions. Equations are expressed as 

(1) 



— = -^V- (pv). 



dt 
d_ 
dt 



(pv) = -V . 



pvv 



B^ 



I 



BB 

47r 



pig^z, 



OB 

-=Vx(vxB), 
dsi 



pi 



Pi^Poi 7 ^ 



(2) 

(3) 
(4) 
(5) 



£l 

PO Cv 

where pi, pi, and si denote the fluctuation of density, pressure, and the specific entropy, v, 
B, 5, and Cv denote fluid velocity, magnetic field, gravitational acceleration, and the heat 
capacity at constant volume. The divergen ce-free condition, i.e. V ■ B = 0, is maintained 
numerically using the method introduced in 



Dedner et al 



(|2002l ). The medium is assumed 
to be an inviscid perfect gas with a specific heat ratio 7 — 5/3. The background values 
of density (po)i pressure (po), and specific entropy (sq) a-re assumed to be in an adiabatic 
stratification with constant gravitational acceleration as 

1/V,d 



Po = Ph 



(6) 



2 



Po = Ph 
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1/V,d-1 



1 I Po 

So = Cv log 

Po 



Hh = 



Ph 



(7) 

(8) 
(9) 



where Pb, Ph, and H^, are the background pressure, density, and pressure scale height at the 
bottom boundary. Vad — (7~ l)/7 denotes the adiabatic temperature gradient. We use the 
relation p — pg + pi, p = pq + pi, and s = sq + si . The redu ced speed of s ound t echnique 

|2012bl ). is used 



Hotta et al 



(hereafter RSST), which is introduced in iRempell ()2005() and 
in this study, and ^ denotes the ratio between the original and reduced speed of sound. 
The investigation of the validity of RSST for the flux rising is given in Hotta et al. (in 
preparation). In the initial condition, the toroidal magnetic flux tube (By) is set near the 
bottom boundary as 

I where r < Id. 

(10) 



By - 



Btexp ( - — 







where r < 2d, 
where r > 2d, 



{X - Xc)^ + (Z - Zc)^, 



Bt 



Snpb 



(11) 
(12) 



/3t + l' 

where /3t = 1.6 x 10^ is the plasma beta of the toroidal field, (xc, Zc) = (l-5i?b, 0.15-ffb), and 
d = O.OSHb are the location and the radius of the flux tube. Two cases with or without 
poloidal magnetic flux sheet are given in A poloidal magnetic flux sheet {Bx) is set 
above the toroidal field as 

B, = Bp exp - ( ] , (13) 



2d + dp, 



/_87rpb_ 
/3p + l' 



(14) 
(15) 



where (3p = 10^" is the plasma beta for the poloidal field, and Zp and dp = O.OSHp are the 
location and thickness of the poloidal layer. The vertical magnetic field is set as B^ — 0. 
Figures[T^ and b show pictures of the initial condition, where the fluctuation of the pressure 
is set as pi — —B'^/{8tt), and si = 0. The radius and strength of the toroidal magnetic 
field is d = 1.8 X 10® cm and _B ~ 1 x 10^ G if we adopt the parameters at the base of 
the solar convection zone; i.e., 6 x 10^ cm, p^ — 0.2 g cm^'^, and the speed of sound 

at the bottom Cb = 2.2 x 10^ cm s~^. These values are plausible at the solar convection 
zone, sinc e the "explosion" p rocess can amplify the magnetic field up to 1 x 10'"' G at 



the base ( Hotta et al 



2012ah . and the magnetic field with this strength can reproduce 



the observational r esult of the tilt angle of sunspot pairs at the surface, i.e. Joy's law. 



(Weber et al 



20111) . Owing to a lack of information on the strength and the distribution 



of poloidal field from observation, we assumed the thickness and the strength of the poloidal 
magnetic field layer as d = 4.8 x 10® cm and Bp ^ 400 G. 

The simulation domain extends as (0,0) < {x/Hb, z/Hh) < (3,1.5). The periodic 
boundary condition is used for all variables at a; = 0, 3i?b- At the top and bottom 
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boundaries, the stress-free and nonpenetrating boundary condition is adopted. A sym- 
metric boundary condition for the magnetic field and density, i.e. dBx/dz ~ dBy/dz — 
dBz/dz = and dpi/dz — 0, is used and entropy is fixed as si = at the top and bottom 
boundaries. The adaptive mesh refinement (hereafter AMR) technique is also used, since 
in this study the magnetic flux is highly l ocalized. A self-sim ilar block-structured grid is 



adopted as the grid of the AMR hierarchy (jMatsumotd . 120071) . In this study, the finest grid 



is equal to the grid using 1024 x 512 in uniform structure. The fourth-order, s pace-centered 



differ ence and fourth-order Runge-Kutta for the time integration are used (jVogler et al 



20051) . The details of our AMR code will be given in our forthcoming paper. 



3. Result and discussion 

Figure [2] shows the result of the case without the mean poloidal field (/3p = oo) aX, t = 
lO^iJb/cb- The magnetic flux tube quickly splits and ceases to ri se. This result is almost the 



same as previous studies ljSchussleiill979tlLongcoDe et al.L 119961 ). Figure[3](red lines) shows 



the temporal evolution of the center of the splitted flux tube {xc, z^) and the circulation F. 
We define the center of a splitted flux tube (xc^c) as the point where the magnetic field 
has maximum value in the region x > 1.5i?b- The circulation F is defined as 

F = / / (V X v)ydxdz, (16) 

where dc — 0.06-ffb- The lifting force caused by circulation and horizontal motion (-Flif) 
and the buoyancy (-Fbuo) can be expressed as 

i^iif = -pVxT, (17) 
Fhuo = "5 / / pidxdz, (18) 

Jxc—dc J Zc — dc 

where Vx is the bulk horizontal velocity of the fiux tube. In the final stage, the lifting 
force (directed downward) and the buoyancy (directed upward) are almost balanced in the 
vertical direction. Since the flux tube rises vertically in the first stage, it is accelerated 
horizontally by the horizontal lifting force pVzT, where Vz is the vertical bulk velocity of 
the flux tube. By this driven horizontal motion, a downward lifting force is induced and 
must be balanced with the buoyancy force at a certain point. 

Figures and d show the distribution of angular velocity and angular momentum. We 
define the center of the rotating motion at the point where vorticity has maximum value 
(see the asterisk in Figl^b). The result shows that the angular velocity increases radially, 
while the angular momentum decreases. We discuss the consequence of such distribution 
in §4. 

Figures [1] and [4] show the result of the case with the mean poloidal field (/3p = 10^°)^. In 
the first stage, similar to the case without the poloidal field, the fiux tube splits and moves 
horizontally with large vorticity (Figs. |4^ and b). Then this vorticity gradually decreases 
(Figs. IH;, d, e, and f). When the flux tube loses almost all circulation the magnetic flux 
begins to rise again (Figs. [T^, and h, |4^, and h). Figure [3] (black line) also shows these 
processes, i.e., re-arisal occurs when circulation decreases around t — 5000i?b/cb- Figure 



^ See also http://www-space.eps.s.u-tokyo.ac.jp/p^7ehotta/movie/twist.html 
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[Th shows that after this process the magnetic flux gets twisted. Although it is not shown, 
these twisted magnetic flux tubes rise to the upper boundary without splitting in the final 
stage of this simulation. We note that, although the fiux tubes have vortices of different 
signs, the signs of the generated twist on two fiux tubes are the same since the directions of 
the poloidal field and the toroidal field are the same between two tubes. Figure [S^, b, and 
c show the values of (V x B)y, By, and (V x B)y/By, respectively. The global structure of 
the poloidal field is shown clockwise, and the average of the values of the twist is positive 
near the core of the toroidal field in both flux tubes. 

We investigated the detailed process of the losing of circulation in the the case with 
the initial poloidal layer. Figure [S] shows the temporal evolution of the magnetic tension 
perpendicular to the magnetic field line {T±), i.e, 

T^=e^-[{B-V)B], (19) 

where ej^ =B x a;/(|B||a;|) is the unit vector perpendicular to the magnetic field and the 
vorticity and w = V x v. The poloidal field lines are suddenly bent (see in Fig. |6]). This 
type of sudden bending event on the magnetic field lines occurs frequently at the point 
where the poloidal magnetic energy is less than the energy of rotation. Since the magnetic 
field is stretched and amplified during this process, the feedback by the magnetic tension 
efficiently suppresses circulating motion. This can be seen in Fig. |6l The high positive 
value of the magnetic tension indicates the feedback effect. We recall that when the mean 
poloidal field does not exist, the angular velocity increases and the angular momentum 
decreases r adially (Fig. [21). This c onfiguration is similar to that of the magneto-rotational 



instability (jChandrasekhail 119601 ) : the rolled- up poloidal magnetic field can be bent up to 
the interior of the fiux tube. When the poloidal magnetic field has finite disturbance, the 
inner part that is rotating faster loses the angular momentum due to the magnetic tension, 
i.e., the interaction with the outer slow-rotating part. In this case, the medium that loses 
angular momentum moves inward radially, and the bending of the poloidal magnetic field 
increases. Therefore the flux tube obtains the twist, which distributes it over the tube's 
cross section without a localization in the boundary layer. We have determined that when 
the magnetic field is bent, the structure of circulation is similar to that of the first case 
without the poloidal field; i.e., the angular velocity increases and the angular momentum 
decreases. 

Our results seem to suggest that two flux tubes may emerge and form active regions in 
different latitudes simultaneously. This probably does not occur because of th e modulation 



of the tube's emerging motion by the thermal convection. IWeber et al.l (j2011l ) have studied 
the rising of the magnetic fiux through the solar convection zone by using the thin-flux-tube 
approximation with the effect of thermal convection and flnd that the timing of emergence 
differs by more than one month between tubes from different latitudes (see their Fig. 8) 
and that the emerging latitude is also scattered typically over ten degrees (see their Fig. 
11) with the strength of a magnetic field of 4 x 10^ G, which is the strength of our splitted 
magnetic flux on their re-rising. Their results suggest that our split pair tubes also emerge 
with these spatial and temporal scattering offsets. Our splitted flux tubes separated about 
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IHp, w hich corresponds to five degrees in latitude and well below the predicted scatter (10 
degrees, IWeber et alj . I2OIII) 



4. Summary 

Using two-dimensional magnetohydrodynamics calculations, we investigated a twist- 
generation mechanism on a magnetic flu x tube at the base of the solar convection zone. It 
is based on the idea of IChoudhuril ([20031) in which a toroidal magnetic field is wrapped by 
a surrounding mean poloidal field. During the generation of the twist, the flux tube follows 
four phases. (1) The flux tube quickly splits into two parts with vortex motions rolling up 
the poloidal magnetic field. (2) Because the physical mechanism is similar to that of the 
magneto-rotational instability, the rolled-up poloidal field is bent and amplified. (3) The 
magnetic tension of the disturbed poloidal magnetic field reduces the vorticity, and the 
lifting force caused by vortical motion decreases. (4) The flux tube gets twisted and begins 
to rise again without splitting. 



With the use of the reduced RSST ( Hotta et al 



2012b() and AMR, it became possible 



to adopt parameters that are acceptable at the base of the solar convection zone, i.e., the 
radius of flux tube d = 1.8 x 10* cm and the strength of the toroidal field i? ~ 1 x 10^ G. In 
this study, we use 400 G as the strength of the poloidal field. It is expected that a certain 
level of strength of such a poloidal field is necessary for suppressing the circulation motion 
and the ensuing flux tube's rising again. The criterion for rising again is an interesting 
issue. The detailed parameter survey will be given in our forthcoming paper. 

Acknowledgements. The authors wish to thank Tomoaki Matsumoto for helping us include AMR in our 
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Fig. 1. Temporal evolution of contours of toroidal magnetic fleld {log[{By / Bq)^]) in cases 
with a poloidal magnetic field, where i3Q/(87r) = phC^- The white rectangles in panels a, 
c, e, and g indicate the region of the panels b, d, f, and h, respectively. The white lines in 
panels b, d, f, and h show the poloidal magnetic fleld lines. 



6 



H. Hotta and T. Yokoyama: Generation of twist 



(a) Magnetic field (c) Angular velocity 




1.97 2.01 2.05 2.09 1.97 2.01 2.05 2.09 

x/Hb x/Hb 

Fig. 2. The results without poloidal field at t — 10*i7b/cb- (a) The contour of toroidal 
magnetic field {\og[{By/Bo)'^]). (b) The vorticity of y-direction in the white rectangle in 
panel a normalized by Cb/i?b- The point where vorticity has a maximum value is indicated 
with an asterisk. This point is defined as the origin of a circular geometry in order to 
calculate angular velocity and angular momentum, (c) The angular velocity normalized by 
Cb/ffb- (d) The angular momentum normalized by pb^bCb- 
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Fig. 3. Temporal evolution of the center of the flux tube, (a) Xc, (b) Zc, and (c) the 
circulation F. The black (red) line denotes the result with (without) the poloidal magnetic 
field. 
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Fig. 5. Values in t = 2000077b/cb- (a) The contour of (V x B)y normalized with Bo/Hh- 
Positive value shows the clockwise poloidal fields, (b) The contour of Toroidal field By 
normalized by B^. (c) The contour of (V x 'B)y/By normalized by Hy^^. The value is only 
calculated in the place where the strength of toroidal field exceeds 6 x 10^'*i?o. 
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(a) t= 3200.00[H,/cJ (c) t= 3600.00[H,/cJ 




Fig. 6. Temporal evolution of magnetic tension {Tj_: Eg. [T^ normalized by Bq/H\^ in the 
case with a poloidal field. Red circles indicate the location of the bending event. 
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